Critical Behaviour in the Relaminarisation of Localised Turbulence in Pipe Flow 
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The statistics of the relaminarisation of localised turbulence in a pipe are examined by direct 
numerical simulation. As in recent experimental data (Peixinho & Mullin Phys. Rev. Lett. 96, 
094501, 2006), the half life for the decaying turbulence is consistent with the scaling (Re c — Re) -1 , 
indicating a boundary crisis of the localised turbulent state familiar in low- dimensional dynamical 
systems. The crisis Reynolds number, is estimated as Re c = 1870, a value within 7% of the 
experimental value 1750. We argue that the frequently- asked question of which Re and initial 
disturbance are needed to trigger sustained turbulence in a pipe, is really two separate questions: 
the 'local phase space' question (local to the laminar state) of what threshold disturbance at a given 
Re is needed to initially trigger turbulence, followed by the 'global phase space' question of whether 
Re exceeds Re c at which point the turbulent state becomes an attractor. 



Understanding the behaviour of fluid flow through a 
circular straight pipe remains one of the outstanding 
problems of classical physics and has continued to in- 
trigue the physics community for more than 160 years 
[3 13 13 13 Although all evidence indicates that the 
laminar parabolic flow is linearly stable, the flow can be- 
come turbulent even at modest flow rates. The exact 
transition point depends not only on the flow rate (mea- 
sured by the Reynolds number Re = UD/v, where U 
is the axial flow speed, D is the pipe diameter and v 
is the fluid's kinematic viscosity) but also sensitively on 
the shape and amplitude of the disturbance (s) present 
0; Q 0; 0- When it occurs, transition is abrupt 
with the flow immediately becoming temporally and spa- 
tially complex. Given that most industrial pipe flows 
are turbulent and hence more costly to power than if 
laminar, a central issue is to understand the conditions 
which trigger sustained turbulence. The problem is, how- 
ever, severely complicated by the fact that the threshold 
appears very sensitive to the exact form of the distur- 
bance and long turbulent transients can exist close to the 
threshold. Of particular interest is the low-Re situation 
where the transition typically leads to a clearly localised 
turbulent structure called a 'puff' within the laminar flow 
0, 0. A puff has a typical length of about 20.D along 
the pipe (see Fig. and, despite appearing established, 
can relaminarise without warning after travelling many 
hundreds of pipe diameters downstream. 

There have been a number of contributions to this 
problem but so far no consensus on the minimum 
Reynolds number, Re Cl above which turbulence is sus- 
tained. Experimental studies have focussed on plotting 
transition-threshold curves in disturbance amplitude- Re 
space for specific forms of applied perturbation. One 
well-studied perturbation having six-fold rotational sym- 
metry gave rise to a threshold amplitude which scaled like 
Re' 1 above Re = 2000 but diverged at Re c « 1800 
0, i.e. below this value no sustained turbulence could 
be excited however hard the flow was disturbed. Sub- 



sequent experiments studying the statistics of relam- 
inarisations of puffs as Re is reduced have lowered this 
threshold value to Re c = 1750 ± 10, close to a previous 
estimate of 1760 but not to others of 1876 10] and 
~ 2000 The only complementary numerical work 
performed so far has been in a short periodic pipe of 5D 
length [ll| where it was demonstrated that the pipe-long 
turbulent state displays the transient characteristics of a 
chaotic repellor until Re c = 2250 above which it becomes 
a chaotic attractor. Recent experiments using a very long 
pipe [l2| in which the statistics on long transients are 
available, however, suggests that there is no critical be- 
haviour. Rather than the turbulent half life scaling like 
r ~ (Re c — Re) -1 0,[U, it is found to increase expo- 
nentially instead. Interestingly, re-interpretation of the 
5D-pipe data seems to corroborate this alternative ex- 
ponential lifetime behaviour even though the pipe is too 
short to capture a turbulent puff. 

In this Letter, we consider a much longer pipe of length 
16ttD (« 50D) in which turbulent puffs can be repre- 
sented faithfully using direct numerical simulation [l^ 
and examine the statistics of how they relaminarise. We 
find an exponential distribution of lifetimes and the crit- 
ical scaling law r ~ (Re c — Re) -1 , with a constant of 
proportionality and an estimate of Re c = 1870 both in 
good agreement with experimental data [3 . Surprisingly, 
given its long history, this represents the first time that a 
quantitative connection between theory and experiment 
has been established in the pipe flow problem. 

The Navier-Stokes equations for an incompressible 
Newtonian fluid, 



d t u - 



Vu + Vp = vV 2 u, 



0, 



(1) 



in a straight pipe with circular cross-section and for con- 
stant mass-flux, were solved numerically in cylindrical 
coordinates (r, z) using a mixed pseudospectral-finite 
difference formulation [2jj. The code was found to 
accurately reproduce linear stability results for Hagen- 
Poiseuille flow, instabilities of nonlinear travelling wave 



FIG. 1: Numerical 'puff' at Re = 1900. (r, z)-section of (V x u) z . Only 20D shown of 50D computational domain. 
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FIG. 2: Probability of relaminarisation after time T at 
Re = 1740 is the same for increased resolutions ('dt' data with 
timestep halved, 'dr' data with 60 radial points, 'dz' data with 
axial resolution of ±576), pipe length ('L' is 100D data), and 
different disturbance ('F' is data obtained with the puff gen- 
erated by an initial period of body forcing - the data is shifted 
in the last case to account for a longer transient period). All 
data sets effectively overlay the default data 'WK'. Inset, nu- 
merical puff spectrum at Re — 1900, A n — maxfc m |A n fe m | 2 , 
index n of Chebyshev transformed radial modes, k, m axial 
and azimuthal Fourier modes respectively and similarly for 
Ak and A m . 



solutions and the statistical properties of turbulent pipe 
flow 14] (as well as being cross- validated with another 
code A resolution of 40 radial points was adopted 

with grid points concentrated at the boundary, Fourier 
modes were kept up to ±24 in 0, and to ±384 in z for a 
periodic pipe of length L = IOttD. This ensured spec- 
tral drop off of 6 orders of magnitude in the power of 
the coefficients when representing a puff velocity field at 
Re = 1900: see inset in Fig. The timestep was dy- 
namically controlled using information from a predictor- 
corrector method and was typically around 0.006 D/U. 
The initial conditions for the calculations were randomly- 
selected velocity snapshots taken from a long puff simu- 
lation performed at Re = 1900. A body forcing applied 
over 10 D of the pipe and for a time 10 D/U was used 
to generate an 'equilibrium' puff which remained stable 
in length and form for a time period of over 2000 D/U 
(see Fig. [Q. At a chosen Re < 1900 a series of at least 
40 and up to 60 independent simulations were performed 
each initiated with a different puff snapshot to generate 
a data set of relaminarisation times. The signature of 
the relaminarisation was a clear and sudden transition 



FIG. 3: The probability of turbulent lifetime > T, P(T), for 
several Re in a periodic pipe of length 16ttD. 



to exponential decay of the energy. The criterion for re- 
laminarisation was taken to be such that the energy of the 
axially-dependent modes was less than 5 x 10~ 4 pU 2 D 3 , 
below which all solutions were well within the decaying 
regime. The range of measured Re c discussed above in- 
dicates sensitivity to noise. Robustness of the relaminar- 
isation statistics was verified by comparing the half-lives 
of data sets obtained by varying different computational 
parameters of the simulation (see Fig. EJ). All modifica- 
tions produce half-life values within the 95% confidence 
interval about the default half-life prediction. 

Decay probabilities for a range of Reynolds numbers 
are shown in Fig. [3| over an observation window of 
1000 D/U. The linear drop-off of the probability on the 
log-plot strongly suggests the exponential distribution, 
P(T) ~ exp(-Tln2/r), where r = r(Re) is the half-life 
of a puff. The median of (T — to) was used as an estimator 
for r. Inspection of the data by varying the cut-off time 
to revealed the effects of an initial transient period in 
the first few data points. This was minimised by select- 
ing to to exclude the first 5-10% of the data (determined 
by looking for the least sensitivity in the half-life predic- 
tion). The results plotted in Fig. B] are consistent with 
the relation r = a(Re c — Re) -1 where a is 2.4 x 10 -4 
compared to 2.8 x 10 -4 obtained in |8( and there is a 
shift of 7% in Re c up to 1870 in the numerical data. 
Also shown is the reinterpreted numerical data for the 
5D pipe and the recent half-life results from the long 
pipe experiments which indicate that 1/r varies ex- 
ponentially with Re rather than linearly. Although the 
data from |l2| is for longer times, there is sufficient over- 
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FIG. 4: The reciprocal of the puff half-life r plotted against 
Re. Data plotted: 'WK'- 50 D data (each data point is the 
result of 40-60 simulations); 'PM' - experimental data from 
0; 'FE' - reinterpreted 5D data 11]; 'H' - experimental data 
from [12. Inset, log-plot of 1/r vs Re. 



lap to suggest that the data from 8] and our results are 
not consistent with being the earlier linear-looking part 
of this exponential. Rather, the results indicate qualita- 
tively different behaviour |28j . 

The exponential probability distribution P(T) found 
here in Fig. |3| implies that puff relaminarisation is a mem- 
oryless process - the probability that the puff will decay 
in a given interval of time is proportional to the length 
of the period but independent of previous events. This 
feature has been found previously in turbulent relami- 
narisation experiments in pipe flow 0,^3 as we ^ 
as in plane Couette flow PJI^PJl an d numerical calcu- 
lations using models of this together with other linearly 
stable shear flows [HI, [2(1 • Faisst and Eckhardt [ll| in- 
terpret this result as indicating that the transient tur- 
bulent state for Re < Re c represents a chaotic repellor 
in phase space. Our results indicate that this conclusion 
carries over to a localised turbulent puff in a long pipe. 
The building blocks for such a repellor are saddle points 
and families of these in the form of travelling waves with 
discrete rotational sym metries are now known to exist 
down to Re = 1251 [2l|, 0, 0. Tentative experimen- 
tal evidence for their relevance to puffs has already been 
found [3 an d corroborating numerical evidence is now 
emerging [l^. The entanglement of all the stable and 
unstable manifolds associated with these saddles at some 
higher Re presumably gives rise to sufficiently compli- 
cated phase dynamics to appear as a turbulent puff in 
real space. That this phase space structure is initially 
'leaky' ultimately allowing escape (relaminarisation) is 
perhaps unsurprising but what is less clear is how it sud- 
denly becomes an attr actor at Re c . The clean scaling of 
the transient decay half life, r ~ (Re c — Re)~ x strongly 
suggests a boundary crisis [25] while the precise value of 
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FIG. 5: Trace of perturbation energy versus additional pres- 
sure fraction required to maintain fixed mass flux, 1 + (3 = 
< d z p > /d z piam (the origin represents laminar flow), for 
the three cases of a sustained puff at Re = 1900 (solid), 
met ast able puff at Re — 1860 with sudden relaminarisation 
(dotted) and the immediate decay of a perturbation (dashed) . 
The inset shows that the energy trace for the metastable puff 
is similar to the sustained puff before it laminarises. 

the critical exponent hints at a simple dynamical systems 
explanation. One, of course, cannot rule out the possi- 
bility that the region never becomes an attractor with 
the exit probability becoming; extremely small but stay- 
ing finite as Re increases [l^. Or, in fact, that there 
are a number of 'leaks' which one by one seal up giving 
a half-life behaviour which varies over a number of dis- 
crete time scales. Also, at some point, the effect of noise 
must surely become significant over long times. However, 
the fact that the numerical simulations and the experi- 
mental results 8] are quantitatively consistent despite 
being subject to different types of errors/disturbances 
indicates that noise is not important over timescales of 
0(1000 D/U) for the levels maintained here and in the 
experiments. 

The simulations confirm that the puff characteristics 
are continuous as Re crosses Re c and that a puff corre- 
sponds to a part of phase space disjoint from the laminar 
state (see Fig. [5] and inset). This observation naturally 
divides the usual question as to how to trigger sustained 
turbulence in pipe flow into two separate issues. Firstly, 
what disturbance at a given Re is needed to trigger tur- 
bulence initially — i.e. what initial conditions will cause 
the flow to leave the neighbourhood of the laminar state 
to reach the puff region of phase space. And secondly, 
what Re is needed so that, for a flow already in the tur- 
bulent region, the flow never leaves — i.e. the puff has 
become an attractor. The implications of this realisa- 
tion are that experimental curves in [5| and [3 showing 
a threshold curve on a disturbance amplitude- Re plot 
must, in fact, be two curves as shown in Fig. El Figure[S| 
shows how initially a threshold amplitude of disturbance 
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FIG. 6: Sketch of the two (independent) thresholds associ- 
ated with transition: one is amplitude-dependent (and highly 
form-dependent) indicating when a turbulent episode is trig- 
gered, and the other is a global i?e-dependent threshold indi- 
cating when the turbulence will be sustained. 



is required to push the solution away from the laminar 
state and into the turbulent region. Once here, the exit 
from the metastable state is sudden and unrelated to the 
entry as relaminarisation is a memoryless process. 

To summarise, numerical simulations described in this 
Letter have clarified the existence of two independent 
thresholds for sustained turbulence. Results probing the 
relaminarisation threshold closely match a recent exper- 
imental investigation 8]. For timescales extending these 
experiments — t < 1000 D/U — we confirm the presence 
of an exponential distribution for the probability of puff 
relaminarisation and corroborate critical-type behaviour 
in which the puff half-life diverges as (Re c — Re) -1 . Good 
quantitative agreement between the experimentally and 
theoretically-estimated value of Re c (less than 7% differ- 
ence) is a rare triumph in this famous canonical problem. 
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